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Abstract. The EGRET instrument has measured detailed 
photon spectra between 30 MeV and 10 GeV, which are rep- 
resented by means of single power-law fits for sources in 
the 3EG catalog. However, various sources show indications 
of spectral cutoffs at GeV energies, which are poorly repre- 
sented by such simple fits. In the case of well exposed or 
bright EGRET sources, a description of spectral cutoffs with 
more complicated functional forms appears to be applicable. 
An application for such multicomponent fits should be seen 
in extrapolations beyond the energies accessible to EGRET, 
i.e. for detectability studies of low-threshold Imaging Atmo- 
spheric Cherenkov Telescopes (lACTs). In cases of unidenti- 
fied gamma-ray sources positionally coincident with Super- 
nova remnants, the spectral shape beyond power-law extrap- 
olations might observationally explain why such prominent 
SNRs hke 7 Cygni, IC 443 and CTAl have not been detected 
during several observation campaigns performed by lACTs 
over the last years. 



1 Introduction 

Spectra of gamma-ray sources in the Third EGRET catalog 
(Hartman et al. 1999) are consistently represented by spec- 
tral indices from single power law fits. In cases of well mea- 
sured source spectra (i.e. PSR J0633H-1746 or PSR B0833- 
45), there are noticeable deviations from single power-law 
model fits reported (Fierro et al., 1997). Here we apply higher- 
order spectral fits as introduced by Bertsch et al. (2000) to 
the dataset of "steady" unidentified gamma-ray sources used 
by Gehrels et al. (2000) for a gamma-ray source popula- 
tion study. This subset of the 3EG catalog sources has been 
used in order to discriminate different source classes among 
the population of unidentified gamma-ray sources, among 
other things by means of a comparison in the (simple power 
law) spectral index. It remains to be investigated, if dif- 
ferences in the average spectral index is a vahd separation 



criteria in terms of the given uncertainties in the spectral fit 
as well as in respect of more appropriate spectral descrip- 
tion using higher order functional forms. Additionally, satel- 
lite based gamma-ray measurements still do not connect to 
the higher neighboring energies, currently accessible using 
ground-based Cherenkov telescopes. Therefore extrapola- 
tions need to be justified, either by appropriate source mod- 
eling and/or incorporating of additional observational infor- 
mation, i.e. a quantitative indication of spectral breaks or 
cut-offs. Here we will especially address cases of unidenti- 
fied gamma-ray sources positionally coincident with SNRs, 
as recently recompiled by Romero et al. (1999) and investi- 
gated at energies above 300 GeV by Buckley et al. (1998). 

2 Analysis 

The "steady" unidentified sources from the 3EG have been 
fitted with three functional forms, the standard single power 
law, two matching power laws, and a power law with expo- 
nential cut-off: 
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Eo has been set to the value determined by the EGRET 
spectral fitting routine as used for the spectral fit given in the 
3EG catalog. The location of the break energy has been set to 
1000 MeV to keep the number of parameters at a minimum. 
With each fit, a reduced was obtained, followed by an 
F-Test to decide if there is statistical justification for the tran- 
sition from a single power law fit to a higher order functional 
form. In the F-Test, a value of p < 0.05 is generally taken as 
a measure where the higher order fit is warranted. Summed 
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data sets for EGRET observations during CGRO observation 
cycles 1 to 4 were taken to maximize the available statistics. 
Spectral studies of EGRET-detected blazars have been pre- 
viously carried out by Mukherjee et al. (1997), Pohl et al. 
(1997) and Lin et al. (1999), but never beyond single power 
law representations. Because spectral hardening has been ob- 
served in gamma-ray blazars during outburst (Sreekumar et 
al., 2001), an overall spectral comparison might not be sim- 
ilarly appropriate for AGN. They will be investigated later 
for different spectral representation than single power law, 
whenever possible on shorter timescales than the four years 
of CGRO observations used here. 



3 Results 

As already reported for the gamma-ray sources close to the 
Galactic plane in Bertsch et al. (2000), the majority of sources 
appears to be best represented by single power law fits. This 
trend has been confirmed for almost all unidentified sources 
at higher Galactic latitudes. However, only for significantly 
detected gamma-ray sources the question of departure from 
a single power law spectrum could be addressed with any 
chance of a statistically meaningful result. Improvements 
from non-Unear fitting are indicated by values lying right 
of the line in Fig.l and Fig. 2. In Fig. 3 the choice of prefer- 
ence of either the two-power law fit or the power law with 
exponential cut-off is indicated by deviations from the line. 
Generally, a more complex spectra could only be established 
if sufficient counts have been observed, rather continuously 
distributed over the energy range accessible for EGRET, and 
obviously with bin entries above 1 GeV. 
However, the F- Tests do not justify the statistical signifi- 
cance of each individual improvement in x^- As indicated 
in Fig. 4 (lines representing the F-Test probability p = 0.05), 
the majority of sources with better representation by non- 
linear fitting functions are insensitive to the form describ- 
ing an spectral cuf-off within the currently available data 
(bottom-left quadrant in Fig.4). A clear indication of a sta- 
tistically significant preference of either the two power law 
fit or the power law with exponential cutoff is only given 
for few individual gamma-ray sources: preference of expo- 
nential cut-off are indicated in bottom-left quadrant, prefer- 
ence of two power law representation in top-left quadrant 
of Fig.4. Although many gamma-ray sources could not be 
compared at all to higher order functional fits due to insuf- 
ficient statistics above the chosen break or cut-off energy, a 
tendency is indicated and will be probed by GLAST: The 
better a gamma-ray source has been observed the more sig- 
nificant deviations from the picture of single power fit repre- 
sentations are needed. EGRET seems to occupy a position 
where instrumental limitations generally do not satisfy more 
complex interpretations of measured source spectra, however 
well-measured individual gamma-ray sources already point- 
ing towards structure in the gamma-ray source spectra unable 
to be adequately represented by a single power-law between 
30 MeV and 10 GeV. 
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Fig. 1. Comparison of the quality of a single power law fit vs. two 
power law fit by means of reduced 



Figure 5-8 give examples of individual unidentified sources 
with different functional representations (Eq. 1-3). All these 
gamma-ray sources are positionally coincident with SNRs, 
three of them suspected of resembling the properties of neu- 
tron stars in gamma-rays (Brazier et al. (1996), Brazier et al. 
(1998), Cheng & Zhang (1998)). These candidates are espe- 
cially interesting because they were studied by ground-based 
Cherenkov telescopes, however yielding only upper Umits so 
far (Buckley et al., 1998). For interpretation of these results 
the model of Drury et al. (1994) has been used to extrapolate 
from the spectrum measured by EGRET. Here we suggest to 
extrapolate from EGRET observational results using better 
spectral representations (Petry & Reimer, 2000). This might 
be an alternative or complementing way, especially under the 
premise that low-threshold Cherenkov telescopes will narrow 
the energetic gap to existing EGRET measurements. 
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Fig. 2. Comparison of tfie quality of a single power law fit vs. power 
law fit with exponential cut-off by means of reduced 



Fig. 3. Comparison of the quality of a two power law fit vs. power 
law fit with exponential cut-off by means of reduced 
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Fig. 4. F-Test probabilities for the significance of the improve- 
ment obtained by using higher order functional form for fitting the 
EGRET spectra. Probabilities less than 0.05 indicate that the more 
complex fitting function is indeed statistically relevant. 
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Fig. 5. Photon energy spectrum of 3EG JOOlO+7309, positionally 
coincident with SNR Gl 19.5+10.2 (CTA 1), shown with three dif- 
ferent representations of a spectral fit. 
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Fig. 7. Photon energy spectrum of 3EG J1801-2312, positionally 
coincident with SNR G6.4-0.1 (W28), shown with three different 
representations of a spectral fit. 
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Fig. 6. Photon energy spectrum of 3EG J0617+2238, positionally 
coincident with SNR G189.1+3.0 (IC443), shown with three differ- 
ent representations of a spectral fit. 



Fig. 8. Photon energy spectrum of 3EG J2020+4017, positionally 
coincident with SNR G78.2+2.1 (7 Cygni), shown with three dif- 
ferent representations of a spectral fit. 



